Josephson Effect Without Superconductivity: Realization in Quantum Hall Bilayers 
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We show that a quantum HaU bilayer with the total filling v = 1 should exhibit a dynamical regime 
similar to the flux-flow in large Josephson junctions. This analogy may explain a conspicuous peak 
in the interlayer tunneling conductance [Phys. Rev. Lett. 84, 5808 (2000)]. The flux-flow is likely 
to be spatiotemporally chaotic at low bias voltage, which will manifest itself through broad-band 
noise. The peak position can be controlled by an in-plane magnetic field. 



Motivation. — In the classic realization of the Joseph- 
son effect, the dynamical variable is the difference be- 
tween phases of pairing wave functions in two supercon- 
ductors. The coupling of this variable to electromag- 
netism is dictated by general principles of gauge invari- 
ance. When pairs tunnel from one superconductor to the 
other, there is a coherent transport of the charge con- 
jugate to this phase difference, and the Josephson effect 
follows. One is led to exactly the same equations, except 
with 2e — > e, when one has a condensate of excitons - 
bound states of particles and holes - extending across a 
junction. This is because an exciton, despite its overall 
neutrality, does couple to the gauge potential difference 
on the two sides of the junction, as a unit charge e. In 
the exciton condensate the coherent charge transfer is 
actually simpler than in a superconductor: it requires a 
single tunneling act instead of two ||l|. Notwithstanding 
the hidden connotations, the excitonic condensate can be 
the ground state of the junction [|| |^. One might be op- 
timistic that in suitable materials exciton Josephson ef- 
fects should be common and could exist at much higher 
temperatures than conventional superconductivity, since 
the fundamental particle-hole interaction is attractive. 

In this Letter we argue that the first observation 
of such effects has already occurred albeit in a low- 
temperature device: a bilayer quantum Hall system with 
total filling factor v = I. For small interlayer separations 
this system spontaneously develops an order parameter 
A = {'4'iip2) oc e*"^, where ipj is the electron creation op- 
erator in layer i [Q . When expressed in terms of fermionic 
operators x\ = i^h creating an electron in layer 1, and 
= ip2, creating a hole in layer 2, A takes the familiar 

BCS form (x|x|) ^-nd enables one to identify (p with the 
phase of the condensate of bound electron-hole pairs. We 
briefly mention that another fruitful analogy is obtained 
by treating the layer index as a pseudospin degree of 
freedom [Q. Along that route, nonzero A signifies pseu- 
dospin ferromagnetism, while (j) determines the direction 
of the in-plane component of the "magnetization" . 

The Josephson effect in the excitonic superfluid was 
originally discussed in Ref. || and more recently in Ref. |^. 
Its appearance in the quantum Hall bilayers was also pre- 
dicted |0, but not immediately accepted Yet a 



complete formal analogy to the Josephson junction does 
hold in the low-energy limit where obeys a damped- 
driven sine-Gordon equation of motion, precisely as the 
phase variable in a large-area Josephson junction. The 
same equation governs many other physical systems, e.g., 
sliding charge-density waves and growing crystals, which 
implies broad relevance of the issues discussed below. 

The aforementioned likely candidate for the Josephson 
phenomenon in a bilayer quantum Hall system is a zero- 
bias peak in the differential tunneling conductance dl/dV 
discovered by Spielman et al. |^ . We interpret this peak 
as an analog of the flux-flow resonance in a Josephson 
junction. The discovery of such a resonance in 1964 by 
Eck et al. was one of the first verifications of the 
Josephson effect in conventional superconductors. The 
experiment of Spielman et al. may well play a similar 
role for the quantum Hall bilayers. 

Two theoretical papers motivated by this experiment 
have appeared Our independent analysis has 

some overlap and some differences with those works. We 
start with a perturbative calculation of I{V). One of 
its predictions is the dependence of the peak position on 
the external in-plane magnetic field. In principle, that 
can be used to map out the dispersion relation of the 
collective mode Q analogous to Josephson plasma os- 
cillations. However, we point out that the perturbation 
theory solution is unstable for low damping characteristic 
of low-temperature quantum Hall systems, which leads to 
a spatiotemporally chaotic fiux-fiow and large tunneling 
current fluctuations. Constructing the theory for such a 
complicated dynamical regime remains a challenge. For 
its numerical study in a discretized formulation (Frenkel- 
Kontorova model), see Ref. |l^ 

Basic equations. — The general form of the equation 
of motion for <f> can be established from considerations of 
symmetry and gauge invariance. To clarify its physical 
meaning we choose another approach. We start with the 
charge conservation equation edtu + d^jk +Jtun = 0. The 
local excess n of the exciton density causes charge imbal- 
ance between the layers. We can express it in terms of 
the local chemical potential difference fj, and the capaci- 
tance c per unit area: n = c/i/e^. The Josephson relation 
entails dtfJ, = hdfcf). The in-plane exciton current j con- 
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sists of the supercurrent jg = {eps/K)dx4> and the con- 
verted quasiparticle current jqp — —dxfJ'/epxx, where ps 
is the condensate phase stiffness and p^x is the resistivity 
of uncondensed quasiparticles. The Josephson tunnehng 
current is jtun = (e/?i)noAsas sin(0 — Qxi), where uq is 
the average electron density per layer, Agas is the tunnel- 
ing strength, Q = eBd/hc is the wavevector imposed by 
the in-plane component B of the external magnetic field 
along the a;2-axis. (Its presence is necessitated by gauge 
invariance). d is the interlayer separation. Assembling 
all the terms, we obtain the equation for (cf. Ref. Q): 

{df - pdtdl -dl)(b + sin((/) + (bo-Qxi)= 0. (1) 

Here we expressed distances in units of "Josephson pen- 
etration length" Xj = (ps/npAsas)^^^ and frequencies 
in units of "Josephson plasma frequency" ujj — v/Xj. 
The velocity parameter v — {e/h){ps/cy^^ ~ OAe^/nh 
is the one that enters the dispersion relation ||] 



V q of the plasmon and parameter /3 = 



{1/A'k){?ujJj / ps){pxxe'^ /h) controls dissipation Fi- 
nally, we included a random phase shift 00 to represent 
disorder (see below). Using the data reported in Ref. ||] 
and theoretical results reviewed in Ref. ^ we arrived at 
rough estimates Aj ~ 5 puv, ujj ~ 10^° s~-^, and /3 ^ 0.01. 
Boundary conditions for Eq. (|l|) depend on the sample 
and measurement geometry. We consider the case where 
each layer is a square with side L, the contact to layer 1 
is along the side xi = —L/2, —L/2 < X2 < L/2, and the 
contact to layer 2 is along the side xi — L/2. This devi- 
ates somewhat from the setup used in Ref. ^ but should 
be inconsequencial as long as the bottleneck for the cur- 
rent flow is the interlayer tunneling not the sheet conduc- 
tivity. This is indeed the case experimentally because the 
inequality V ^ Ipxx is satisfied {V is the voltage differ- 
ence between the contacts). Under such conditions it is 
also permissible to choose (/>i — (pR — Vt, and use 



I{t) = / dx2[dMt) - ML{t)] 



(2) 



to calculate the tunneling current. Here (f>L,R = 4>i^i — 
±L/2,X2) and our units of voltage and electric current 
are Vg = hujj/e and Ig = eps/h, respectively. 

Before we proceed to the calculations let us explain the 
origin the the term "flux-flow." Let us consider the case 
(j)o = dt(f> = V = and focus on the limit Q ^ 1 where 
the stable (ground-state) solution of Eq. (P corresponds 
to an almost uniform phase distribution 1. There 

is an equivalent alternative formulation in terms of a 
shifted phase = ip — Qxi in which Q disappears from 
the argument of the sine but re-emerges in the bound- 
ary conditions for diO = dip — Q. In this formulation 
the ground-state 6 varies rapidly and almost linearly in 
space, which can be described as the 27r/(3-periodic lat- 
tice of 27r-solitons. A nonequilibrium state for V ^ Q 
where the phase increases uniformly with the rate V can 



then be visualized as a uniform sliding of the soliton lat- 
tice — whence the term flux- flow. 

Perturbation theory. — The perturbative solution of 
Eq. (U) for 00 = is readily done in terms of the Green's 
function 



G{uj, k) = [lj"^ + ia{k)uj - 



Y{k) = Pk^. (3) 



The dc current / is then obtained by averaging Eq. (g) 
over time. The full expression is somewhat cumbersome 
but for L in the range 1 <C L ^ /^^ it reduces to [M 



/ = 



aiQ)V 



2 (V^2 _ Q2)2 +Q,2(Q)y2 



(4) 



We estimate L ~ 50 and / 10 ^ for the conditions of 
Ref. |, so L is in the required range. From Eq. M) we see 
that the resonance arises when the velocity V/Q of the 
soliton train matches the plasmon velocity v (v = 1 in 
the adopted dimensionless units). Under this condition 
the power dissipation in the system reaches a maximum. 
More generally, the resonance condition is eV = hujo{Q) 
where loq^Q) is the Agas — > limit of the plasmon disper- 
sion relation (in which the plasmon becomes the Gold- 
stone mode). The last conclusion was reached indepen- 
dently in Refs. |ll| and ^ 

Let us turn to the disordered case, (po ^ 0. From 
Eq. (|^) and the identity sin z = Im e'^ , one can see that 
the random phase factor e**^" plays the same role for the 
disordered system as the periodic phase factor e''^^'^ for 
the clean system. The Fourier expansion of e'^^° can be 
thought of as a set of "diffraction gratings" with different 
wavevectors fc, resonating whenever V — k. As a result, 
the resonance is broadened in proportion to the breadth 
of the Fourier spectrum of the following correlator: 



a(x) = (e-'"^o(x)gz0o(o)^^ 
and the appropriate generalization of Eq. (^) is 

d^k 



I = Im 

2 



:a(k-Q)G{V,k). 



(5) 



(6) 



To evaluate the integral we assume that the random 
phase field <po is mainly due to static randomly posi- 
tioned and randomly oriented bound vortex pairs in the 
exciton condensate. In the context of the pseudospin 
ferromagnet analogy mentioned earlier such vortices 
were recognized as merons, the topological defects of the 
0(3) nonlinear cr-model, and a remarkable fact was es- 
tablished: each meron possesses an overall electric charge 
e/2 concentrated near its core. The competition of the 
Coulomb repulsion between the like electric charges and 
attraction between oppositely charged vortices selects the 
optimal size a of meron pairs; a remains finite below the 
roughening transition ||l^ temperature Tr = Sirps/kB. 
From our point of view, in the context of the Joseph- 
son junction analogy, the meron pairs correspond to mis- 
aligned Abrikosov vortices trapped in the junction. The 



I-V characteristic of such junctions was investigated by 
Fistul and Giuliani These authors also attempted 

to go beyond the perturbation theory but that part of 
their calculation is suspect for reasons that will become 
clear shortly. 

Within the chosen model of disorder we obtain 



CT(fc) Aa^{kay-^, 
A ~ 27r7, 
A ~ 27rln(fca)"\ 



7 = Trn^a 
7 = 2, 



(7) 
(8) 
(9) 



where rim is average density of the meron pairs. Substi- 
tuting the above formulas into Eq. (||) yields for Q = 
and K < ^ < /?" 



-1. 



I = I,{V,/Vf-'', I, = L^a{V,)V, 



(10) 



Here K = [A(V;)a^/8/3]i/(^-'^) is determined from the 
condition that the variance of 0, which is given by 



(2^ 



a(k-Q)|G(T/,fc)|^ 



(11) 



becomes of the order of one. AtV < V, the premise of the 
perturbation theory (0 ^ 1) is violated and formula ( p^ ) 
is invalid. However, y = 14 is not the true boundary 
of the perturbative domain. Below we will show that 
the perturbation theory actually breaks down in a wider 
range of voltages, V < K*, where K* ^ K for small p. 

In the following we focus on the case 7 « 2 (not an 
unreasonable value [pT|) where the Fourier spectrum of 
e"^" is extremely broad, see Eq. (|^). The flux-flow reso- 
nance is equally broad so that it shows up as a plateau 
in I{V), see Fig. ^ With the chosen value of 7 we get 
the expression (physical units temporarily restored) 



(12) 



for the height of the plateau. Similar expressions were 
also derived in Rcfs. |l^ and |l^. 

Even though the response at low V is beyond the reach 
of the perturbation theory, the "Josephson critical cur- 
rent," which is a zero-bias parameter, can nevertheless 
be calculated. The collective-pinning theory of Vinokur 
and Koshelev [nsl yields (for 7 = 2) 



I, = 2ieps/h){L^ay\j)ln'{Xj/a) 



(13) 



Reminiscent of the inequality between the coefficients of 
static and kinetic friction, Ic exceeds /*. Note that Ic 
is essentially a static quantity. It need not coincide with 
the y — > limit of the dynamical response formulas akin 
to Eq. (^. The overall dependence of / on F is likely to 
appear the way shown in Fig. [l[ 

Choosing 7 « 2 enables us to reproduce the plateau 
feature in the experimental data |^ (the peak in dl/dV 
corresponds to a plateau in I vs V), but the theoretical 
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FIG. 1: I-V characteristic given by the perturbation theory 
with 7 ~ 2 (solid lines). The dashed line is a conjecture. 
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FIG. 2: Kinematic diagram for the parametric resonance. 



estimate of is off by two orders of magnitude. How- 
ever, this estimate is highly unreliable in view of large 
uncertainties in the parameters of the model and of their 
strong renormalization by thermal fluctuations, known 
from the theory of the roughening transition . 

Instability of a tachyonic flux-flow. — It is known, al- 
though not widely appreciated, that the perturbation 
theory solution for the case (j)o — - in the form of 
a moving soliton lattice - can become unstable when 
the lattice velocity exceeds the velocity of Josephson 
plasmon, V/Q > 1. The instability can be under- 
stood as a parametric resonance driven by the exter- 
nal source of frequency V and wavevector Q. Indeed, 
if I^ > Q, then kinematics allows a simultaneous exci- 
tation of two counter-propagating plasma waves, with 
wavevectors fc± = (Q ± V)/2, see Fig. ||. 

For large Q the instability appears when (cf. Ref. Eol) 



(14) 



and so the perturbation theory solution on the "tachy- 
onic" side V > Q is stable only if V > 14* = {-i/PV^^ 
{Q <^ K* is assumed). This condition is much more re- 
strictive than the naive V — Q ^ derived from the 
criterion \(f>\ <^ 1. 

Very close to the threshold 14*, the parametric res- 
onance described above produces a small modulation of 
the uniform soliton train |l^ , but as V moves closer to Q 
other instable modes proliferate and the system dynam- 
ics quickly enters the regime of spatiotemporal chaos. In 
Fig. H we show the comparison between the results of 
the naive perturbation theory and of our numerical sim- 
ulations of a one-dimensional system. As expected, they 
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FIG. 3: Top: I-V characteristic. Dots are from numerics, 
solid line is from a full perturbation theory formula for a finite- 
length system. It contains both the main resonance of Eq. (^) 
and smaller "Fiske resonances." The scatter of the datapoints 
at V > 2.5 is due to statistical noise. The parameters of 
simulation are L = 20, Q — 3, P — 0.01, and 00 ~ 0. Bottom: 
An example of the frequency spectrum of I{t) in the above 
simulation. Note the narrow peak at u equal to the applied 
voltage V — 4 and the strong broad-band noise at smaller 

Lj's. 



agree when either V > 14* or Q — V > Q^^, but at 
intermediate V^s they differ substantially and a strong 
broad-band noise in / sets in. 

Let us now show that a similar instability must occur in 
the disordered case as well. To estimate the correspond- 
ing we approximate the unstable modes by wavepack- 
ets of plasma waves with a Gaussian spread of wavevec- 
tors, ?A±(k) = (27rAfc2)-icxp[-|k-k±|V2(Afc)2]. If Afc 
is smaller than a(fc±), then G(fc,k) is approximately the 
same for dominant k's within the wavepackets, and so the 
left-hand side of the instability criterion (|lj) need not be 
modified. The right-hand side however is changed from 
1/4 to {l/2n)a{V - Q){Aky and the threshold voltage 
becomes P^^^^a^^^^ . Thus, for small /? the lowest 

voltage at which the perturbation theory is expected to 
apply is y** S> 14. Our conjecture on the behavior of 
I in the interval K < F < (depicted in Fig. |^) is 
based on a notion that chaos leads to a higher effective 
dissipation, sec Fig. ^ (top) and Refs. |Q and The 
suppression of / below is an educated guess motivated 
by the work of Fistul and Giuliani [|l6| . 

The probable observation of Josephson effects without 
bulk superconductivity in a bilayer quantum Hall ferro- 



magnet paves the way for exploring a vast variety of other 
Josephson phenomena in this system. It should also in- 
spire attempts to realize them in other systems. The 
complexity and appeal of the emerging theoretical issues 
(not fully appreciated in the current literature), as well as 
their recurrence in other contexts warrant further study. 
From this perspective the lack of quantitative agreement 
between the perturbation theory and experiment is stim- 
ulating. 
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